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ABSTRACT 

QUaD is a bolometric CMB polarimeter sited at the South Pole, operating at frequencies of 100 and 150 GHz. 
In this paper we report preliminary results from the first season of operation (austral winter 2005). All six 
CMB power spectra are presented derived as cross spectra between the 100 and 150 GHz maps using 67 days 
of observation in a low foreground region of approximately 60 square degrees. This data is a small fraction 
of the data acquired to date. The measured spectra are consistent with the ACDM cosmological model. We 
perform jackknife tests which indicate that the observed signal has negligible contamination from instrumental 
systematics. In addition by using a frequency jackknife we find no evidence for foreground contamination. 
Subject headings: CMB, anisotropy, polarization, cosmology 
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1. INTRODUCTION 

The CMB is expected to be polarized at the ~ 10% level 
due to Thomson scattering by free electrons of the local 
quadrupole in the CMB radiation field at the time of last scat- 
tering. The resulting polarization signal can be decomposed 
into two independent modes. At the time of last scattering, 
even parity E-modes are generated by both scalar and ten- 
sor (gravitational wave) metric perturbations while odd-parity 
B-modes are generated only by gravitational waves. A sec- 
ondary source of B-mode polarization comes from the weak 
gravitational lensing effect of intervening large scale struc- 
ture, which converts E-modes into B-modes on small scales. 

The first detection of the E-mode polarization signal 
was made by the 30 GH z radio interferometer, DASI, in 
2002 dKovac et all 12002)) . Since then, in a ddition to a 
furthe r measurement by the DASI experiment dLeitch et al.l 
2005), E-mode measurements have been made with the 
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CBI dReadhead et alJl2004h. CA PMAP dBarkats et al.ll2 005). 
BOO MERanG dMontrov et all |2006l) . WMAP dPageetalJ 
2006) and MAXIPOL (|WjieLalj20O6) experiments. 

High precision measurements of the E-mode signal rep- 
resent a non-trivial test of the standard cosmological model 
since the polarization of the CMB probes the velocity field 
at the time of last scattering, as opposed to the density field 
probed by CMB temperature measurements. In addition, ac- 
curate measurements of E-mode polarization can be useful 
for constraining certain cosmological parameters which are 
fairly insensitive to the CMB temperature field (e.g. isocur- 
vature modes in the early Universe). Such a high resolution 
measurement of the E-mode polarization signal is the pri- 
mary science goal of the QUaD 1 experiment. In addition to 
the polarization, QUaD will, with further analysis, also pro- 
vide interesting results on the CMB temperature field on small 
scales. 

In this paper we present preliminary power spectra mea- 
sured from QUaD's first season of operation. The paper is or- 
ganized as follows. In Section|2] we summarize the QUaD in- 
strument and describe the observation strategy and low-level 
data reduction. Section[3]describes our map-making and sim- 
ulation procedure. Our power spectrum estimation method 
is described in Section [4] and the power spectrum results are 
presented in Section [5] along with results from a number of 
jackknife tests. In Section [6] we estimate cosmological pa- 
rameters using our spectra and our conclusions are presented 
in Section [7] 

Throughout this paper when we refer to "the ACDM 
model" we mean spe cifically the model genera ted by the 
CMBFAST program (IZaldarriaga & Seljakl l2000h using the 
WMAP3 cosmological paramet ers given under the heading 
"Three Year Mean" in table 2 of lSpergel etall d2006l) . This is 
the model used in our simulations and shown in the plots. 



2. INSTRUMENT SUMMARY AND OBSERVATIONS 

QUaD is a millimeter-wavelength bolometric polarimeter 
designed for observing the CMB at two frequency bands, 100 

1 QUaD stands for "QUEST and DASI". In turn, QUEST is "Q & U 
Extragalactic Survey Telescope" and DASI stands for "Degree Angular Scale 
Interferometer". The two experiments merged to become QUaD in 2003. 
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and 150 GHz. The experiment is sited at the MAPO obser- 
vatory, approximately 1 km from the geographic South Pole. 
First light was achieved in February 2005, and science ob- 
servations began in May 2005. The telescope is a 2.6 m on- 
axis Cassegrain with nominal beam sizes of 6.3 (100 GHz) 
and 4.2 (150 GHz) arcminutes. The tower, ground shield and 
altitude-azimuth mount of the DASI experiment are re-used 
for QUaD, the ground shield being extended to accommo- 
date the larger telescope structure. The mount has a third axis 
which rotates about the optical symmetry axis (termed "deck" 
rotation). This is a very useful feature for a polarimeter, as it 
allows the entire telescope to be rotated to an arbitrary angle 
with respect to the sky. 

The QUaD receiver comprises two anti-reflection coated 
cryogenic re-imaging lenses and a focal plane array of 3 1 pix- 
els, each composed of a corrugated feed horn and two orthog- 
onal P olarization-Sensitive Bolometers (PSBs; iJones et all 
(2003)). The PSB pairs are oriented on the focal plane in 
two groups with bolometer sensitivity angles separated by 
45 degrees. This redundancy in detector orientation allows 
one to construct maps of the sky in Stokes Q & U with ob- 
servations at a single deck angle if so desired. Each pixel is 
single-frequency and the pixels are divided between the two 
observing bands with 12 at 100 GHz and 19 at 150 GHz. The 
PSBs are similar to those flown on the successful B2K ex- 
periment dMasi et al.ll2006l) . A complete description of the 
receiver along with details of the optical testing and charac- 
terization will be provided in Hinderks et al. (in prep). 

The first season of QUaD observations were completed in 
October 2005 and consisted of ~ 1 00 days of CMB runs, in ad- 
dition to special runs for pointing model determination, beam 
mapping, and detector time constant measurements. The 
CMB runs consist of scanning the telescope back and forth 
by 7.5 degrees in azimuth, in a series of 30 second constant- 
elevation "half-scans". The telescope is then stepped by 0.02 
degrees in elevation and the process repeated to build up a 
raster map. Since the telescope is sited close to the Earth's 
axis of rotation, azimuth and elevation closely approximate 
to right ascension and declination. In the first season we 
have mapped a 60 square degree pa tch in the low-foreground 
B2K deep region dMasi et al.l l2006). For our chosen scanning 
speed and observing declination I ~ 2000/, where / is the 
frequency in Hz at which multipole number I appears in the 
time ordered data. The timeconstants of most (80%) of our de- 
tectors are less than 30 ms with the slowest two ~ 100 ms cor- 
responding to half power roll-offs at £ ~ 10000 and I ~ 3000 
respectively. 

To permit the removal of ground contamination, the scan- 
ning strategy employs a "lead-trail" scheme whereby each 
hour of observations is split equally between a "lead" field 
(1st half hour) and a "trail" field (2nd half hour), separated by 
0.5 hours in right ascension. The lead and trail field obser- 
vations follow exactly the same pattern in telescope azimuth 
and elevation so a constant ground signal can be removed by 
differencing the lead and trail field data. Furthermore, each 
day of observation is split into two 8 hour blocks made over 
different ranges in azimuth with the telescope rotated at dif- 
ferent deck angles. This enables a powerful jackknife test as 
described below. 

Initial processing of the raw time-ordered data (TOD) con- 
sists of glitch removal, deconvolution of the bolometer tem- 
poral response, low pass filtering and down-sampling. The 
relative calibration factor between channels (and within chan- 
nels of a given pair) is derived from frequent short scans in 



elevation (el-nods) which introduce a strong atmospheric gra- 
dient into the TOD. This relative calibration is applied sepa- 
rately within each frequency group. Various quality control 
data cuts are applied at this stage: days with bad weather, bad 
pointing, poor focal plane temperature stability, or moon con- 
tamination are discarded. After applying these data cuts, 67 
of the ^100 days of CMB observations remain, and are used 
in the following science analysis. 

3. MAP-MAKING AND SIMULATION PROCESS 

Two analysis pipelines have been constructed which are in- 
dependent in the sense that they share no code, although the 
algorithms are intended to be identical (with some important 
exceptions — see below). For this initial analysis we use data 
which has been point by point lead-trail differenced as de- 
scribed above. There is clear ground pickup in the data, which 
although mostly common mode, has a polarized component. 
This pickup appears to be completely canceled in the field dif- 
ference. We note that many CMB experiments have mitigated 
ground pickup by field differencing, for example DASI. 

Before mapping, a best-fit 3rd order polynomial is sub- 
tracted from each half-scan to remove the bulk of the 1/f 
noise. The signal from each PSB pair is then summed and dif- 
ferenced and co-added into the map according to the telescope 
pointing information. In the co-addition, a weighting is ap- 
plied according to the inverse variance of each 30 second half- 
scan to properly account for differences in sensitivity across 
pairs, and to down weight periods of poorer weather. The pair 
sum data is used to construct total intensity (T) maps and the 
pair difference data is used to construct maps of Stokes Q and 
U by using the known orientation angle of each detector pair 
with respect to the sky. During this process, a ~ 10% cor- 
rection for the non-ideal polarization efficiency is applied to 
the differenced PSB data. These angles and efficiencies are 
measured from special observations of a chopped, polarized, 
thermal source placed externally to the telescope. 

While building the T, Q and U maps, we also construct the 
expected T, Q and U variance maps under the simple assump- 
tion that the noise is uncorrected (white) in the TOD. These 
variance maps are used later to weight the signal maps in the 
power spectrum estimation stage. 

As described in the next section the power spectrum esti- 
mation technique requires signal only, noise only, and signal 
plus noise simulations of the complete experiment. We gen- 
erate these as follows. 

To construct simulated noise timestreams, we Fourier trans- 
form each set of half-scans and take the auto and cross spectra 
between all channel pairs. This allows us to re-generate sim- 
ulated noise timestreams with the same frequency spectra as 
the real data and the proper cross correlations between chan- 
nels. We assume that the instantaneous signal to noise in the 
time-ordered data is sufficiently low that the power spectra 
created in this way accurately estimate the noise. 

To generate simulated signal timestreams we generate real- 
izations of T, Q and U sky maps under the ACDM model us- 
ing the "synfast" generator (part of the HEALPix package 2 ), 
convolve with the instrumental beam (separately for each de- 
tector), and re-sample according to the pointing of the tele- 
scope. Included in this process are a scatter in the polariza- 
tion efficiency factors of the PSB pairs and a scatter in the 
PSB sensitivity angles. The uncertainties used are derived 
from the special observations of a chopped, polarized ther- 

2 See http://healpix.jpl.nasa.gov/index.shtml and Gors ki et al] 42005) 
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mal source mentioned above. A detailed beam model is used 
which is derived from special beam mapping runs on the com- 
pact HII region RCW38, daily scans of each detector across 
this source, and observations of a bright quasar PMN J0538- 
4405 which lies within our field. 

Either of the signal or noise simulated timestream, or their 
sum, is then passed through the standard mapping algorithm 
(complete with polynomial subtraction and variance weight- 
ing) to yield simulated maps. 

To derive the absolute calibration factor of our e xperiment 
we pass the B2K temperature maps ( Ma si et al.l2006l) through 
the simulation process to provide maps which are filtered in 
exactly the same way as the QUaD maps, and are thus di- 
rectly comparable to them. Both sets of maps are then Fourier 
transformed and cross spectra taken between them to deter- 
mine the relative calibration factor. B2K is in turn calibrated 
against WMAP. (WMAP3 lacks sufficient sensitivity within 
our small sky area to allow a direct cross calibration.) We per- 
form this calibration separately for the 100 GHz and 150 GHz 
maps, resulting in 5% absolute calibration uncertainty in units 
of temperature (10% in units of po wer), including b oth QUaD 
and BOOMERanG uncertainties dMasi et al.ll2006l) . To moni- 
tor the absolute stability of the instrument an internal calibra- 
tion source is inserted into the beam frequently during routine 
data taking and these readings show excellent stability over 
the entire season of observations (few percent for any given 
channel and 0.5% for the gain ratio of a PSB pair). 

For visual presentation only, we divide by the variance map, 
transform to E and B, and filter to the angular scales where 
signal to noise is highest, to produce Figure Q] 

4. POWER SPECTRUM ESTIMATION METHOD 

To estimate angular power spectra, we employ a Monte 
Carlo (MC) based analysis. This method requires the creation 
of noise only simulated power spectra to correct the measured 
power spectra for noise bias, and signal only power spectra to 
allow the suppression of power by filtering to be corrected. In 
addition signal plus noise spectra are required to provide the 
final covariance matrix of the bandpower measurements. 

Before measuring power spectra from the real and simu- 
lated maps, we apply an inverse variance mask to the maps 
based on the expected spatial distribution of the noise (using 
the variance maps mentioned above), and additionally mask a 
small number (5) of point sources which are apparent in our T 
maps, and confirmed by external catalogs. We note that none 
of these sources is detected at high significance in the Q and 
U maps. 

At this point the two pipelin es diverge — one f ollows the 
standard MAST ER technique oflffivon et al]d2002l) . extended 
to polarization dBrown et alJl2005l) . and works explicitly on 
the curved sky (pseudo-Q), while the other makes the flat sky 
approximation and uses two dimensional Fourier transforms 
to derive power spectra. 

The first pipeline measures raw pseudo-Q power spectra 
from the maps using a modified version of the "anafast" pro- 
gram included in the HEALPix package. Estimates of the 
CMB power spectra are then reconstructed as bandpowers 
(Pi) from the pseudo-Q measurements using 
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, = ^M^^P^(C,-(N,) MC ). 



(1) 



Here, Pbg is a binning operator in ^-space, Ci are the raw 
pseudo-Q spectra measured from the real data and (N^ ) M c are 
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FIG. 1. — QUaD first season field differenced polarization maps decom- 
posed into E and B-modes, and filtered to include only the angular scale range 
200 < I < 1000. The top row shows the result for signal (non-jackknife) 
maps, while the bottom row is for the "deck" jackknife (see text). This plot 
shows 150 GHz maps and the color scale is ±30 (i¥L- in all cases. 



the average pseudo spectra measured from the noi se only sim- 
ulation s. Mbb' is the binned coupling matrix of iBrown et al.1 
d2005l) which corrects for mode-mixing due to the survey ge- 
ometry. M.bf also contains the correction for the effects of 
both the TOD polynomial filtering and the telescope beam 
width. These corrections are derived from the set of sig- 
nal only simulations. Finally, the covariance matrix of the 
resulting bandpowers is found from the scatter among the 
power spectra measured from simulations containing signal 
and noise. 

The second pipeline takes the two dimensional Fourier 
transform of the masked maps, converts the Q and U Fourier 
modes into E and B and calculates bandpowers as the mean of 
the product of the modes within each annular bin. The product 
can be taken as the auto spectrum of a given map, or as a cross 
spectrum between two maps (which need not be at the same 
frequency). This is done for the real data maps and for each 
simulation realization. The data spectra then have the mean of 
the corresponding set of noise only simulations subtracted to 
noise correct them. Filter/beam suppression factors are calcu- 
lated as the ratio of the mean of the signal only simulations to 
the input power spectra, and the data spectra are corrected by 
dividing out these suppression factors. Finally the bandpower 
covariance matrix is estimated from the scatter of the signal 
plus noise simulations in the same way as for pipeline one. 

The first pipeline explicitly corrects for mixing between the 
EE and BB spectra due to the sky cut using the ma- 
trix. The second pipeline does not make such a correction, 
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but the level of mixing into the BB spectra in the signal only 
simulations is found to be negligible compared to the current 
instrumental sensitivity (< 0.2 /xK 2 ). In addition note that the 
simulations indicate that inter- and intra-spectra mode mixing 
due to the half-scan polynomial filtering, and lack of cross 
linking in the map, are absolutely irrelevant for the current 
analysis. 

For either pipeline we can take spectra internally within the 
sets of 100 and 150 GHz maps, or cross spectra between the 
two frequencies. In this paper, we present only the frequency 
cross spectra. The single frequency spectra will be presented 
in a future paper. 

5. RESULTS AND JACKKNIFE TESTS 

Figure [2] shows the frequency cross spectra measured 
from the first season of QUaD data by the two independent 
pipelines 3 . The error bars are the square root of the diagonal 
elements of the bandpower covariance matrices which are es- 
timated from the run-to-run scatter among MC simulations as 
described in the previous section. The first pipeline includes a 
mode decoupling step which narrows the £ range to which 
each bandpower responds, while at the same time chang- 
ing the adjacent bandpower correlation coefficients from their 
"natural" value of « +0.2 to w -0.2. This leads to an increase 
in the diagonal of the bandpower covariance matrix which is 
reflected by larger error bars on the plot. However we empha- 
size that the total information content of the bandpowers from 
both pipelines is similar. In Figure [3] our results are shown 
compared to published results from other experiments. 

A powerful test for systematic contamination from ground 
pickup, or another source, is the so-called jackknife test. In 
this paper we use a map based jackknife, forming separate 
maps from various approximately equally sized data subsets, 
subtracting these and taking the power spectra of the result. 
We also do this for the signal plus noise simulations to esti- 
mate the expected uncertainty of the jackknife spectra. In as 
much as the signal originates on the sky it should exactly can- 
cel under jackknife — depending on the split, and its origin, 
false signal is not likely to do so. Here we use the follow- 
ing data splits each of which will be explained in turn: deck 
split, scan direction split, season split, focal plane split and 
frequency split. 

The so called deck split is possibly the most powerful test. 
As mentioned above each day of observations on our CMB 
field is split into two 8 hour blocks. Because the run starts 
always at the same local sidereal time these blocks occur al- 
ways over the same range of azimuth angle as the telescope 
turns within its ground shield. In addition each block of ob- 
servations is made always at the same rotation angle of the 
telescope with respect to the line of sight, with a 60 degree 
rotation occurring between the two blocks. Hence each given 
bolometer pair scans the sky at a different orientation angle 
within each block. Therefore the deck jackknife polarization 
maps will only cancel if the rotation to the absolute refer- 
ence frame which occurs in the map making step is being per- 
formed correctly. The ground pickup is observed to be very 
complex with certain pair differences showing a detectable 
spike always at a certain azimuth angle and rotation angle of 
the telescope. Hence even if some ground pickup is leaking 
through the field differencing operation, we certainly do not 
expect it to appear identically in the deck split Q and U maps 
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FIG. 2. — Preliminary QUaD first season results derived as cross power 
spectra between the 100 and 150 GHz maps. The two sets of points are the 
results from two completely independent analyzes of the data (see text for 
details) and are displaced by ±10 in multipole number from their nominal 
values for clarity. 



Bandpowers, covariance matrices and bandpower window functions are 
available in numerical form at http://quad.uchicago.edu/quad 
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FIG. 3. — Comparison of preliminary QUaD first season power spectrum 
results to selected other publishe d resu lts to date. The experiments com- 
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and thus to cancel in the deck jackknife. Figure |4] compares 
the signal (non-jackknife) and deck jackknife power spectra 
and we see that the vast majority of the apparent sky signal 
cancels. Note that where the signal spectra are sample vari- 
ance dominated the jackknife spectra error bars are smaller 
since there is no sample variance in a null spectrum. 

For the scan direction jackknife we form separate maps 
from the half-scans in each direction. If deconvolution of the 
detector temporal response is not done correctly then the for- 
ward and backward maps will not match and residuals will 
remain when they are subtracted. In some cases the time con- 
stants of the two halves of a detector pair are not well matched 
and hence poor deconvolution could lead to leakage from T 
into polarization making this a very important test. 

The split season jackknife forms maps from the first half 
and second half of the used days. If for example there were 
a drift in the absolute calibration of the instrument over time 
then cancellation failure would be expected here. 
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FIG. 4. — QUaD signal spectra compared to deck jackknife spectra. (See 
text for details.) 
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TABLE 1 





PTE 


VALUES 


FROM X 2 


TESTS 






Jackknife 


TT 


TE 


EE 


BB 


TB 


EB 


Deck angle 


0.236 


0.208 


0.812 


0.435 


0.274 


0.062 


Scan direction 


0.000 


0.173 


0.304 


0.375 


0.236 


0.223 


Split season 


0.032 


0.814 


0.257 


0.527 


0.904 


0.111 


Focal plane 


0.193 


0.702 


0.079 


0.503 


0.450 


0.225 


Frequency 


0.034 


0.306 


0.610 


0.452 


0.642 


0.135 


Signal" 


0.501 


0.964 


0.415 


0.482 


0.066 


0.809 



a The PTE value for the signal case is calculated against the ACDM 
model. 

The focal plane split forms separate maps using the two ori- 
entation groups of bolometer pairs in the focal plane. Because 
observations are taken at two deck angles it is possible to con- 
struct Q and U maps using each group. 

For the frequency jackknife instead of taking cross spec- 
tra between the 100 GHz and 150 GHz maps we subtract 
them and take the spectra of the frequency difference maps. 
Any admixture of two or more signal components with dif- 
fering spatial distributions and frequency spectral indexes is 
expected to fail this test (for example CMB plus synchrotron 
and/or dust) — this is therefore a stringent test for foreground 
contamination. Figure[5]compares the signal (non-jackknife) 
and frequency jackknife power spectra and we see that to the 
limits of experimental sensitivity the sky pattern is identical 
at the two frequencies. 

Figures|4]and|5]are visually impressive but to quantify how 
well these tests are passed we have calculated \ 2 statistics for 
the comparison of the jackknife power spectra with the null 
model. For some of the jackknifes we do not expect perfect 
cancellation due to the interaction of the polynomial filtering 
by the half-scan and the non perfectly overlapping coverage 
region of the two data subsets. Hence we compare the mea- 
sured x 2 values to the distribution which we measure from the 
set of signal plus noise simulations rather than to a theoretical 
X 2 distribution. 

In Table Q] the probability to exceed (PTE) the x 2 value 
measured from the data is tabulated for each spectrum and 
jackknife test. These numbers are for the first pipeline — the 
second pipeline gives results which are similar. Ideally these 
PTE values should be uniformly distributed from zero to one. 
The only value which shows an obvious problem is in the 
scan direction case — the enormous signal to noise of the TT 
spectrum seen in Figures [4] and [5] makes the TT jackknifes 
sensitive to tiny systematic errors. The jackknife tests of the 
single frequency spectra show some additional problems pass- 
ing these formal \ 2 tests and we are hence choosing not to 
publish these spectra at this time. 

The final row of Table [TJ is not a jackknife — it is a com- 
parison of our measured spectra against the WMAP3 ACDM 
model mentioned above and shown in the plots. The PTE val- 
ues lie within the acceptable range indicating that our mea- 
sured spectra are consistent with this model. 

6. COSMOLOGICAL PARAMETER ESTIMATION 

We have carried out a basic 6 parameter cosmological 
parameter constraint analysis using our polarization power 
spectra only, i.e. we use the TE, EE and BB spectra, but 
not TT. Our methodology uses the Monte-Carlo Markov 
Chain (MCMC) me thod and is based o n that of the WMAP 
team as described in lVerde et al.l d2003l) : we impose the same 
flat priors, the same re-parameterization and the same con- 
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FIG. 5. — QUaD signal spectra compared to frequency jackknife spectra. 
(See text for details.) 
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TABLE 2 
Polarization only cosmological 
parameter constraints using 
preliminary quad first season 
power spectra 



Parameter 



Symbol 



Value 



Baryon density fl/,h 2 

Matter density Q m h 2 
Hubble constant h 
Optical depth r 
Scalar fluctuation amplitude 3 A s 
Scalar fluctuation index" n s 



0260 +0 0061 
u.uzou _ (| mm 

n 1 49 +0 029 
u -^ z -0.031 

0.74 + 016 
156 + ^ 9 
0.78 



0.833 



a The pivot point for A s and n s is k p = 
0.05 Mpc" 1 

vergence/mixing test. In addition we optimize the Markov 
chain step size choice, a fundamental parameter for the cor- 
rect behavior of the algorithm, by calculating the parameter 
covariance ma t rix of a preliminary run of our chains as in 
Tegm arket al.l (|2004). We deal with the b eam and calibration 
uncertainty using a method proposed by iBridle et"ail (120021) 
which accomplishes an effective marginalization over these 
parameters by adding extra terms to the bandpower covari- 
ance matrix. 

The theoretical p ower spectra are calculated using CAMB 
dLewis et al.ll2000l) . and then transformed by means of the ex- 
perimental bandpower window functions to predictions for 
the binned f/,'s. We assume that the likelihood distribution 
for our bandpowers is Gaussian, which our signal plus noise 
simulations indicate is a valid assumption. The bandpower 
covariance matrix for the TVs is assumed to be independent 
of the model, and is estimated from the signal plus noise sim- 
ulations as described above. 

In Table |2] we present the marginalized expectation values 
and 68% confidence limits for each parameter. Four con- 
verged MCMC chains with around 50,000 steps each are 
merged to obtain these results. For all of the parameters 
our 68% confidence limit encloses the WMAP3 expectation 
value. The \ 2 of the model with the parameters listed in the 
table is 44. 1 giving a probability to exceed this by chance of 
0.51 (for the 45 degrees of freedom). 

7. CONCLUSIONS 

We have presented preliminary power spectra measured 
from the first season of observations with QUaD. In this pa- 
per we have presented only the frequency cross spectra taken 



between the 100 and 150 GHz maps. We find that these spec- 
tra are entirely consistent with the ACDM model — the mea- 
sured EE spectrum has a distinctive peak at I ~ 400 exactly as 
expected, the TE spectrum shows the expected correlations, 
and the BB spectrum is consistent with zero. A basic polar- 
ization only parameter constraint analysis yields confidence 
limits which agree with the WMAP3 results, and are as tight 
as those which were derived from CMB temperature spectra 
just a few years ago. 

We have performed jackknife tests by measuring power 
spectra from differenced maps generated under several data 
splits and find that the results are free from significant in- 
strumental systematics. In addition we have presented a fre- 
quency jackknife which indicates that contamination of the 
CMB by astrophysical foregrounds is negligible for the cur- 
rent experimental sensitivity. 

We note that this analysis considers only frequency cross 
spectra and includes only 67 days out of a total of 250 days of 
CMB observing to date. 

The QUaD experiment has begun a third season of obser- 
vations and analysis of the second season data is underway. 
When completed, we fully expect to improve substantially on 
the preliminary results presented here. 
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